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Abstract 

In the left-right twin Higgs model, one of the neutral Higgses is a natural candidate for WIMP 
dark matter. We analyzed the dark matter relic density in this framework and identified regions 
of parameter space that provide the right amount of dark matter. We also studied the dark 
matter in the more general inert Higgs doublet model in which the mass splittings between the 
dark matter and other particles do not follow the relations in the left-right twin Higgs model. 
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I. INTRODUCTION 



The Standard Model (SM) of particle physics has been extremely successful in explaining 
the data from particle physics experiments. It is, however, not satisfactory when it comes 
to cosmology One of the puzzles that remain is the amount of cold dark matter in the 
Universe, which has been determined precisely from WMAP[l|: 

tt DM h 2 = 0.112 ±0.009. (1) 

None of the SM particles can be a good candidate for the dark matter. The existence of 
the dark matter itself points to new physics beyond the SM. Weakly interacting massive 
particles (WIMPs) are promising dark matter candidates among various proposals in the 
literature. Given the natural size of its interactions (weak interaction) and natural value 
of its mass around the TeV scale, the WIMP relic density is around the observed amount 
of dark matter without much fine tuning. 

Recently, the twin Higgs mechanism has been proposed as a solution to the little 
hierarchy problem [2, 0, Higgses emerge as pseudo-Goldstone bosons once a global 
symmetry is spontaneously broken. Gauge and Yukawa interactions that break the global 
symmetry give masses to the Higgses, with the leading order being quadratically divergent. 
Once an additional discrete symmetry (twin symmetry) is imposed, the leading one-loop 
quadratically divergent terms respect the global symmetry. Thus they do not contribute 
to the Higgs masses. The Higgs masses do obtain one-loop logarithmically divergent 
contributions, resulting in masses around the electroweak scale when the cutoff is around 
5-10 TeV. 

The twin Higgs mechanism can be implemented in left-right models with the discrete 
symmetry being identified with left-right symmetry|3|]. In the left-right twin Higgs (LRTH) 
model, several physical Higgs bosons remain after the spontaneous symmetry breaking. In 
particular, there is an SU(2) l Higgs doublet h, which only couples to the gauge boson sector 
due to some additional discrete symmetry. The lightest particle in its neutral components 
is stable, therefore a good candidate for WIMP dark matter. 

The dark matter candidate in the left-right twin Higgs model is similar to the one in 
the inert Higgs doublet model (IHDM)|5]. In the IHDM, one of the Higgs doublets plays 
the role of the SM Higgs doublet, which couples to both the fermion matter sector and 
the gauge boson sector. The other Higgs doublet, however, only couples to the gauge 
boson sector similar to h in the LRTH model. Therefore, it could contain a dark matter 
candidate. 

There are some differences between these two models though, which are relevant for the 
dark matter studies. In the IHDM, the mass for the SM Higgs boson is a free parameter, 
which can be anywhere between 100 GeV up to a very large value. The SM Higgs mass 
in the LRTH model, on the other hand, is more constrained. It depends weakly on other 
model parameters and is typically between 160 — 180 GeV. The Higgs potential in the 
LRTH model is also more constrained compared to that of the IHDM. The dominant 
contributions to the Higgs potential in the LRTH model come from the one-loop Coleman- 
Weinberg (CW) potential. The quartic couplings due to gauge interactions are about the 
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size of g 4 /(167r 2 ), much smaller than typical quartic couplings in the IHDM, in which 
the natural size is 0(1). In the LRTH model, however, there are extra heavy states and 
interactions, which lead to new features in the relic density analysis that are absent in the 
IHDM. 

The relic density analysis of the dark matter in the IHDM has been studied in Ref . 0] . 
Their results agree with our analysis. We have furthermore identified regions of parameter 
space that have previously been overlooked in Ref. [6|. In addition, due to the extra 
particles and interactions that appear in the LRTH model, there are new regions that 
could produce the right amount of relic density, which are otherwise absent in the IHDM. 

This paper is organized as follows. In Sec. [TTl we briefly introduce the LRTH model, 
emphasizing the Higgs sector. In Sec. IHH we identify the dark matter candidate and 
introduce the parameters that are involved in the analysis. In Sec. IIVt we present in 
detail our numerical results for the dark matter relic density. We discuss two dark matter 
mass regions: (A) low mass region, and (B) high mass region, and identify the regions 
of parameter space in which we obtain the right amount of dark matter relic density. In 
Sec. El we conclude. 



II. THE LEFT-RIGHT TWIN HIGGS MODEL 



The LRTH model was first proposed in Ref. [3| and the details of the model as well as 
the Fevnman rules, particle spectrum, and collider phenomenology have been studied in 
Ref. [7|. Here we briefly introduce the model and focus our attention on the Higgs sector. 

In the LRTH model proposed in Ref. |3j, the global symmetry is U(4) x U(4), with 
the diagonal subgroup of SU(2) L x SU(2) fi x U(1)b-l gauged. The twin symmetry which 
is required to control the quadratic divergences of the Higgs mass is identified with the 
left-right symmetry which interchanges L and R. For the gauge couplings #2L and g2R of 
SU(2)i and SU(2)#, the left-right symmetry implies that g2L = 92R = Qi- 

Two Higgs fields, H and H, are introduced and each transforms as (4,1) and (1,4) 
respectively under the global symmetry. They can be written as 

Hl \ £j ( H L 



where Hl r and Hl,r are two component objects which are charged under SU(2)i x 
SU(2) R x U(1) B -l as 

H L and H L : (2, 1, 1), H R and H R : (1, 2, 1). (3) 

H couples to both the gauge boson sector and the fermion sector while H couples to the 
gauge boson sector only (for reasons discussed below). This can be achieved by imposing 
a discrete symmetry under which H is odd while all other particles are even. Each Higgs 
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acquires a non-zero vacuum expectation value (vev) as (/ 3> /) 



/0\ 



< H > 








< H > 








(4) 
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breaking one of the U(4) to U(3), respectively. Each yields seven Nambu-Goldstone bosons 
and one massive radial mode. The Higgs vevs also break SU(2)jj x U(1)b-l down to the 



Below the cutoff scale A, the radial modes are integrated out and the effective theory 
can be described by a nonlinear sigma model of the 14 Goldstone bosons. We follow the 
parametrization of the scalar fields H and H as described in Ref. {zj. After spontaneous 
global symmetry breaking, three Goldstone bosons are eaten by the massive gauge bosons 
Wft and Zh, and become their longitudinal components. After electroweak symmetry 
breaking, three additional Goldstone bosons are eaten by the SM gauge bosons W ± and 
Z. With certain re-parametrizations of the fields, we are left with four Higgses that couple 
to both the fermion sector and the gauge boson sector: one neutral pseudoscalar 0°, a pair 
of charged scalars (j)^, and the SM physical Higgs hsu- In addition, an SU(2) L doublet 
h = (hijh®), which resides in H L , is charged under an additional discrete parity and 
couples to the gauge boson sector only. 

The left-right discrete symmetry ensures that the global symmetry is respected for 
the one-loop quadratic order of the Higgs potential, and so the quadratically divergent 
corrections to the masses of the Goldstone bosons (i.e., the Higgses) are absent. The sub- 
leading contribution is only proportional to In A, for A being the cutoff scale. No severe 
fine tuning is introduced for A of the order of 5—10 TeV. 

Note that in the LRTH model, it is necessary to introduce both the Higgs fields H 
and H, with a hierarchy between their vevs: / ^> /. The reason is that the masses of 
the heavy gauge bosons, which are proportional to gf, need be large enough to avoid 
the constraints from the electroweak precision measurements. Such a large value for / 
does not reintroduce the fine tuning problem for the Higgs mass, since the gauge boson 
contributions to the Higgs mass are suppressed by the smallness of the gauge couplings. 
The Higgs field H, however, cannot couple to the SM fermions, in particular, the top quark 
sector. Otherwise, the heavy top quark (the fermionic top quark partner that is introduced 
in the LRTH model) obtains a much larger mass of the order of yf, which could result 
in a large, dangerous contribution to the Higgs mass from the top sector. To avoid this, 
a parity is introduced in the model under which H is odd while all other fields are even. 
This parity thus forbids renormalizable couplings between H and fermions, especially the 
top quark sector. Therefore, H couples only to the gauge boson sector, while H couples to 
both the gauge sector and the matter fields. The SU(2)^ doublet h, which resides in Hl, 
is also odd under this additional matter parity. The lightest particle in h, typically one of 
the neutral components, is stable, and therefore constitutes a good dark matter candidate. 



SM U(l) 
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III. DARK MATTER CANDIDATE 



The masses for hf and h 2 , which are relevant for the dark matter relic density analysis, 
can be obtained from the one-loop C vVpotential arising from gauge boson contributions. 
At leading order, the Higgs masses are|7|: 



rn 



2 3 rn 2 V[] 
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2n A 2 ^ 2g 2 + g 2 n A 2 ^ 
g 2 2 (ln + 1 + y \ y2 In — + 1 

m w H 2 m z H 



(5) 



<-< + T ^^(f-l^ + ^ + l), (6) 

where mw,w H ,z,z H are the gauge boson masses. The gauge couplings g\ and g 2 are related 
to e and Weinberg angle 9 W by 

V cos 26 ) l() sin ^ 

The cutoff scale A is typically taken to be Airf. Numerically, the masses for hf and h 2 
are around 200 — 700 GeV for typical values of mw H , m z H i n the model. The small mass 
splitting between hf and h® [2nd term in Eq. (jSJ)] is around a few hundred MeV, which is 
due to electromagnetic interactions that are absent for h®. 
In addition, we can introduce a term in the Higgs potential 

V H = {PhIHl, (8) 

which gives additional masses for hf and h® as 

K ~ Am| 2 ~ fi 2 . (9) 

This allows us to vary the mass of the dark matter independently as a free parameter. 
However, such a term breaks the left-right symmetry softly. Therefore, it is natural for it 
not to be much bigger than / (typically between 500 GeV and 1500 GeV). Note that since 
fi 2 could have either sign, the total mass for hf and h® could be either larger or smaller 
than the contributions from the CW potential in Eqs. (jSJ) and (Q. 
The complex scalar h 2 can be written as 

h° 2 = (10) 

where 5 and A are the scalar and pseudoscalar fields, respectively. The cross section 
for scattering between S (or A) with matter via Z exchange (see Fig. [T]) is about eight 
to nine orders of magnitude larger than the limit of spin independent interactions from 
current dark matter direct detection experiments [8|]: agj < 10~ 42 cm 2 at 90% C.L.. Such 
constraints, however, can be avoided if there is a small mass splitting between S and A: 
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FIG. 1: Feynman diagram that contributes to the elastic scattering of the dark matter with 
nuclei. 

82 = — nig > a few hundred MeV, given the typical kinetic energy of the dark matter 
and the momentum transfer between the dark matter and the scattering nuclei in such 
scattering processes. Here and in our numerical analysis below, we assume that < m^, 
therefore S is the dark matter candidate. For m§ > m^, the pseudoscalar A becomes the 
dark matter candidate. However, the numerical results that we present for S also apply 
for the pseudoscalar dark matter. 




H L A 5 H L m\ L H L 

FIG. 2: Contribution to the Higgs mass term m^^^H^ from the A5 quartic term in the Higgs 
potential. 

Such a mass splitting between S and A, can be obtained by introducing a quartic term 
in the Higgs potential: 

V H = -^[{H[H L f + h.c.]. (11) 

The sign of the quartic term is picked such that > wig for positive A5. Once Hl obtains 
a vev (0,v/\/2), it generates a splitting between m| and m\: 

m\-m\ = \ 5 v 2 . (12) 
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Introducing such a term in the Higgs potential breaks the left-right symmetry explicitly. 
Therefore, it could give a potentially dangerous quadratically divergent contribution to 
the Higgs mass. However, quartic terms —(X^/4)h 2 S 2 and (X 5 /4)h 2 A 2 have opposite signs, 
which ensures the cancellation between the quadratically divergent contribution from the 
S loop and that from the A loop. This leaves only logarithmic contributions which are 
safe for A5 ~ 1. Another way to understand this cancellation is that at one- loop level, 
a four point interaction (\ 5 /2)(H\ j H L ) 2 (both H L fields flow out) cannot give rise to the 
operator m 2 HL H\H L (with one H L flows in and another H L flows out) unless there is a 

mass insertion of m 2 HL H\HL to flip the field flow (see Fig. [2]). Dimensional analysis tells us 
that we can at most get m 2 HL log A terms from such a diagram while the A 2 contribution is 
absent. Given that A5 is almost not constrained, we can get a large splitting between rrig 
and via such a A5 quartic interaction. 

Due to the small mass splitting (of about a few hundred MeV) between and rn^, 

introducing a mass splitting between S and A also leads to a mass splitting between 
the charged Higgses and the dark matter candidate S (or A, if rrig > m£). Letting 
Si = m £ — rrig, we have S2 ~ 25i for 8 2 ~ 1 GeV. On the other hand, we could also 
introduce a quartic term: 

V H = \a\H{H l \ 2 . (13) 

If Hl obtains a vev (0,v/y/2), it generates a mass splitting between ra| and m| (in 
addition to the mass splitting introduced by the U(l) em in the CW potential): 

2 2 ^4 2 /1 A\ 

m? — m t = v . 14 

hi hi 2 

In the LRTH model, such a term also breaks the left-right symmetry which could produce a 
dangerous contribution to the Higgs mass if A4 is too large. Unlike A5, which is not strongly 
constrained as discussed above, a simple dimensional analysis shows that |A 4 | < (//A) 2 . 
For the cutoff scale A taken to be 4nf, | A4 1 < l/(167r 2 ). Therefore, in the LRTH model, 
the corresponding mass splitting between hf and h 2 is at most about a few GeV, typically 
much smaller than 5 2 - Including A4 and A5 terms, and the U(l) em contribution in the CW 
potential, we get 



ml — m| = — v 2 — —v 2 + small CW contribution. (15) 
h! s 2 2 



In the LRTH model, when A5 could be large while A4 is constrained to be small, the 
approximate relation S2 ~ 2Si holds for not so small values of 82- In the general IHDM, 
however, such an approximate relation could be violated. In our analysis below, we study 
both the most general case that 81 and 82 are two independent parameters, as in the general 
IHDM, as well as the case when the relation 82 ~ 28 1 holds, as in the LRTH model. 

IV. DARK MATTER ANALYSIS 

We analyzed the relic density using the program MicrOMEGAs[9]. This program solves 
the Boltzmann equation numerically, using the program CalcHEP[10|] to calculate all the 
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relevant cross sections. The CalcHEP model files for the LRTH model can be found at 



Ref . [ljj . When the mass splittings between the dark matter candidate and other particles 
are small, coannihilation effects are also included. We have identified two mass regions of 
m§ that could provide the amount of dark matter relic density that is consistent with the 
WMAP result at the 3 a level: (A) low mass region where < 100 GeV, and (B) high 
mass region where 400 GeV < < a few TeV. 



A. Low Mass Region: m$ < 100 GeV 



8^100 GeV 8 2 =100 GeV 




(a) 



FIG. 3: Plot of the relic density as a function of nig for the low mass region. The left plot (a) 
fixes Sx at 100 GeV while S 2 is set to be 5 GeV (solid), 10 GeV (dashed), 15 GeV (dotted), and 
100 GeV (dash-dotted). The right plot (b) fixes 5 2 at 100 GeV while Si is set to be 5 GeV (solid), 
10 GeV (dashed), 13 GeV (dotted), and 100 GeV (dash-dotted). The horizontal band shows the 
3 a region of the WMAP result: tth 2 = 0.112 db 0.027. 



Let's first study the most general case when 5i and 5 2 are two independent parameters. 
In Fig. [3X we show the relic density dependence on the mass of the dark matter m§ for (a) 
#1 = 100 GeV (when coannihilation between S and hf can be ignored), and various values 
for S 2 , and (b) #2=100 GeV (when coannihilation between S and A can be ignored), and 
various values for 5i. 

For 5i = 100 GeV and 5 2 = 10 GeV [dashed curve in Fig. [3] (a)], coannihilation between 
S and A needs to be taken into account due to the relatively small mass splitting 5 2 . There 
are three mass windows where Qh 2 falls within the WMAP bounds. The first two mass 
windows near m§ ~ 40 GeV receive dominant contributions from the coannihilation of 
SA — > qq/ll via Z exchange, as shown in Fig. HI The contribution from Zh exchange is 
suppressed by the heavy Zh mass, while the contribution from 0° exchange is suppressed 
due to the small SA(p° coupling. The dip between these two mass windows is due to the Z 
pole: m§ + = m z . We denote these mass windows as "low mass Z pole region". 
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FIG. 4: Feynman diagrams that contribute to the coannihilation processes: SA — ► qq/ll. 




FIG. 5: Feynman diagrams that contribute to the annihilation processes SS — > VFIF/ZZ. 



Once mj increases above the Z pole region, the coannihilation cross section decreases, 
leading to a relic density Qh 2 above the WMAP region. However, when is close to 
the W boson mass, a new annihilation channel SS — > WW opens, which brings the relic 
density down to the WMAP region [the third mass window: nig ~ 75 GeV, of the dashed 
curve in Fig. [3] (a)]. We denote this mass window as "low mass WW bulk region". The 
corresponding Feynman diagrams are shown in Fig.0 When rug further increases, channel 
SS — > ZZ also becomes accessible, which further increases the annihilation cross section, 
reducing the relic density. When 2m § = rrif lSM ~ 170 GeV (as in the LRTH model), the 
dominant contribution comes from the exchange with an Hsm P°l e an d the relic density 
reaches its minimum. The slight rise in Qh 2 around 90 GeV is due to the destructive 
interference between the SS — > WW/ZZ Feynman diagrams. 

For smaller value of 82 , coannihilation effects become more important. The solid curve 
in Fig. [3] (a) shows the Qh 2 dependence on for 62 = 5 GeV. The contribution from the 
coannihilation cross section becomes so large that only part of the low mass Z pole region 
survives. 

As 82 is increased, coannihilations becomes less effective. The dotted curve in Fig. [3] (a) 
shows the relic density for 82 = 15 GeV, when coannihilations can be ignored. The low 
mass Z pole region no longer exists for such large values of 82- The relevant annihilation 
process for the low mass region (before annihilation into WW/ZZ opens up) is SS — > qq/ll 
via h S M exchange (Fig. EJ). The cross section is small due to the suppressed coupling of 
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FIG. 6: Feynman diagrams that contribute to the annihilation processes SS — ► qq/ll. 



SShsM and the small Yukawa couplings of hsMqq and hsM^, which leads to a relic density 
larger than the WMAP preferred value. The low mass WW bulk region, however, still 
remains once m§ is close to the W boson mass. 

As 62 further increases, process SS — > qq/ll becomes more and more important due 
to the enhanced coupling SShsM, which is proportional to 62- The dash-dotted curve 
in Fig. [3] (a) shows the relic density for 62 = 100 GeV, when the relic density falls into 
the WMAP window for a smaller value of m^: m§ ~ 50 GeV. The dominating process 
for this region is SS — > bb. We denote this mass window as "low mass bb bulk region", 
distinguishing from the low mass WW bulk region, which appears for smaller values of 
62- Typically, for <5 2 < 50 GeV, SS — > WW is the dominating annihilation process, which 
prefers a larger m§ around 70 GeV for the WMAP window. For 5 2 >50 GeV, SS — > bb 
starts to dominate, when the WMAP mass window shifts to smaller values of m§. As we 
will show later, the switching between these two regions occurs at lower 62 once Si becomes 
smaller. There is, however, no clear separation of the low mass bb bulk region and the low 
mass WW region since in general both processes contribute to dark matter annihilations. 

Part of the small m§ region, however, has already been excluded when we impose Z 
width constraints. For m§ + < mz, process Z — ► S + A contributes to the Z decay 
width Yz- On the other hand, Y z has been measured very precisely at the LEP[l3|. which 
agrees very well with the SM prediction. Therefore, there is little room for any new Z 
decay process. The region that corresponds to 2m $ + S 2 < m z is excluded once we impose 
the Z width constraints. Therefore, part of the low mass Z pole region is ruled out by the 
Z decay width consideration. 

Fig. [3] (b) shows the relic density as a function of ma for various splittings of Si, with <5 2 
fixed at 100 GeV. The coannihilation effects between S and A can now be safely ignored, 
while coannihilation effects between S and have to be taken into account when Si is 
small enough. We obtain results similar to those presented in Fig. [3] (a). The low mass 
pole region appears for a lower m^ region since it corresponds to coannihilation process 
Shf — > qq' via W exchange instead of Z exchange. We indicate this as "low mass W pole 
region". For larger values of Si ^ 15 GeV, the coannihilation between S and h^ 1 can be 
neglected, and the low mass W pole region disappears. 
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Due to the large value of S 2 , the low mass bb bulk region appears for rrig ~ 50 GeV. 
The low mass WW bulk region does not appear for 62 = 100 GeV since SS — > bb already 
saturates the annihilation cross section at smaller rag. 

The charged Higgs has been searched for at the LEP and the Tevatron 13, 13]. A 
lower mass bound of 74 — 79 GeV at 95% C.L. is obtained at the LEP Il3j considering 
H + — > cs, H + — > t + u, or H + — > W*A. A more recent search at CDF[14j studied the 
charged Higgs produced in the top quark decay t — > H + b, with H + further decaying into 
a pair of quarks, leptons, or W + (p. All of those searches, however, rely on the couplings 
of the charged Higgs with fermions, which is absent for hf in both the LRTH model and 
the IHDM. Therefore, the experimental limit on the charged Higgs mass does not apply. 
On the other hand, Z — > hfh± is allowed once the charged Higgs mass is below 
which leads to the deviation of Z decay width V z from the SM predictions. Therefore, 
= nig + 5i < m z /2 is excluded by the consideration of rz[l2|, and part of the low 
mass W pole region in Fig. [3] (b) is not allowed given this constraint. 




FIG. 7: Plot of the relic density as a function of rug for the low mass region, with 5± and 82 
chosen as: (5 1 ,S 2 ) = (10 GeV, 10 GeV) (solid curve), (10 GeV, 20 GeV) (dashed curve), (10 
GeV, 80 GeV) (dotted curve), and (80 GeV, 10 GeV) (dash-dotted curve). The horizontal band 
shows the WMAP 3 a region. 



Fig. [7] shows the dependence of Qh 2 as a function of trig for four different sets of (Si, 62). 
When 5i and S 2 are both small, (Si, S 2 ) = (10 GeV, 10 GeV), as indicated by the solid curve, 
the coannihilations between S, A and hf should all be taken into account. Therefore, there 
are two dips for around 35 — 40 GeV, which correspond to the W pole and the Z pole, 
respectively. There is a relatively wide region of that falls into the WMAP 3 a region 
for 50 GeV ^ < 70 GeV, which is a combination of the pole region and the WW bulk 
region. For (^i,^) = (10 GeV, 20 GeV) (dashed curve), only the low mass W pole region 
appears since 62 is large enough that the coannihilation between S and A is ineffective. 
However, 62 is still relatively small such that the bulk region is dominantly SS —>■ WW. 
When 5 2 is set to be 80 GeV (dotted line), the bulk region is dominantly SS — > bb due to 
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the relatively large value of <5 2 . For (^i,^) = (80 GeV,10 GeV) (dash-dotted curve), only 
the low mass Z pole region (due to the coannihilation of S and A) and low mass WW bulk 
region appear. 




FIG. 8: Contour plots of the relic density in the nrig vs. 82 plane [plot (a)] and rrig vs. 5\ plane 
[plot (b)] for the low mass region. 5\ (82) is chosen to be 10 GeV (solid lines) and 100 GeV 
(dashed lines) for plot (a) [plot (b)]. The shaded region is excluded by experimental constraints 
from Tz- The region enclosed by two contour lines corresponds to the WMAP 3 a region. 



Fig. M (a) shows the contour plot of Qh 2 in the m§ vs. 62 plane for two values of 81. 
The region enclosed by two contour lines corresponds to the WMAP 3 o region for the 
dark matter relic density. The shaded region is excluded by the LEP constraints on the Z 
decay width (We also applied Y z constraints for a given Si.), which eliminates part of the 
low mass W or Z pole region for small m§. 

For Si = 100 GeV (regions enclosed by the dashed contour lines), the coannihilation 
between S and hf can be ignored. For 7 GeV< 8 2 ^$13 GeV, there are two possible mass 
windows of m§ for a given 62 , corresponding to the low mass Z pole region and WW bulk 
region as discussed earlier for Fig. [3] (a). For 82 ^13 GeV, there is only one mass window of 
rrig for a given 62- This corresponds to the case when the coannihilation effect is negligible 
and the dominating contribution comes from annihilation between dark matter themselves. 
For 82 ~ 50 GeV, annihilation into WW dominates, which leads to mj around 70 GeV 
for the WMAP preferred region. For larger 82, annihilation into bb final states starts to 
contribute. The curve bends to the left, leading to a WMAP mass window at a smaller 
value of m,g. 

For 81— 10 GeV [solid curves in Fig. [S] (a)], the vertical region around rrig = 40 GeV 
corresponds to the the low mass W pole region. The bulk region shifted to the left 
comparing to 8\ = 100 GeV case. For small 82 around 10 GeV, the coannihilations between 
S, A, and hf are effective. There is a relatively large window of rug between 40 to 75 GeV 
that could give the right amount of relic density. 
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Fig. El (b) shows the contour plot of Qh 2 in the m s vs. Si plane for S 2 =100 GeV (dashed 
curves) and 10 GeV (solid curves). For larger values of Si > 15 GeV, the coanihilation 
between 5* and hf can be ignored. The m s window is independent of S± since either the 
Z pole region (for smaller 8 2 ) or bulk region (WW or bb) is independent of Si. For large 
5 2 =100 GeV (region enclosed by the dashed curves), annihilation of SS — > bb results in a 
bulk mass window of m s ~ 50 GeV. When 5 2 gets smaller, the bulk mass window shifts 
to larger m s as SS —>■ WW becomes more and more important. For 5 2 = 10 GeV (regions 
enclosed by the solid curves), coannihilation of S and A needs to be taken into account, 
which leads to the appearance of the low mass Z pole region for m s around 45 GeV. 
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FIG. 9: Qh 2 vs. m s plot (a) and relic density contour plot (b) for 5 2 = 2S±. Si is taken to be 5 
GeV(solid), 10 GeV(dashed), and 20 GeV(dotted) in plot (a). The band in plot (a) and region 
enclosed by two contour lines in plot (b) are the WMAP 3 a regions. 



In the LRTH model, the mass splittings Si j2 follow the approximate relation S 2 = 2Si for 
not so small values of 5 2 . The corresponding relic density plots are shown in Fig. [9] (a) for 
the m s dependence. For small Si j2 , the coannihilation between S, A, and hf leads to the 
low mass pole region. The low mass bulk region is due to SS — > WW. The corresponding 
WMAP 3 a region contour plot of the dark matter relic density is shown in Fig. [9] (b). 



B. High Mass Region: 400 GeV< m s < 1 TeV 

When the dark matter mass m s > 400 GeV, there are regions of parameter space that 
could also provide the right amount of dark matter relic density. Fig. [10] (a) shows the 
relic density dependence on m s for various values of Si. There are typically two mass 
windows. A "high mass bulk region I" appears near 450 GeV < nig < 750 GeV, with 
the dominating contributions coming from SS — > WW/ZZ. A "high mass pole region" 
appears for nig around 810 GeV, where the coannihilation of 5 and A via the exchange of 
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FIG. 10: Plot (a) shows the relic density Q/i 2 vs. rug for 82 = 1 GeV. Plot (b) shows the contour 
plot of the relic density in the mg vs. 62 plane. <5i is chosen to be 1 GeV (solid curve), 5 GeV 
(dashed curve), and 10 GeV (dash-dotted curve) for both plots. / is chosen to be 600 GeV. The 
band in plot (a) and region enclosed by two contour lines in plot (b) are the WMAP 3 a regions. 



a heavy gauge boson Zh dominates (see Fig. H] for the corresponding Feynman diagram) . 
When 8\ gets larger, the SS annihilation cross section grows, which leads to a suppressed 
Qh 2 below the WMAP window. However, at larger nig > 850 GeV [see dash-dotted curve 
in Fig. [10] (a)], a new bulk mass window, denoted as "high mass bulk region II", might 
appear. The numerical results of the dark matter relic density for fixed 8% and varying 62 
are very similar to those presented in Fig. [TU](a). The relic density decreases for increasing 
82 as well. 

Fig. [10] (b) shows the WMAP 3 a window in the m§ vs. 82 plane for various values of 
8\. The pole region around m§ ~ 800 GeV is independent of 8\ t 2- For small values of 8\ t 2, 
the high mass bulk region I always appears. The high mass bulk region II, m§ > 850 GeV, 
on the other hand, only appears when either 8\ or 8 2 is around 7—15 GeV. For 5 12 > 15 
GeV, the annihilation cross sections become so big that the relic density always falls below 
the WMAP preferred region. 

Note that the pole region only appears in the LRTH model, where the existence of the 
Zh pole dominates the coannihilation SA — > qq/ll at m§ ~ mz H /2. In the general IHDM, 
the pole region disappears. However, the general results that we obtained above for the 
high mass bulk regions still apply. 

Fig. [Til shows the WMAP 3 a region in the trig vs. 82 plane when we impose the LRTH 
relation 82 = 28±. A relatively large region of parameter space could provide the right 
amount of dark matter relic density. 

The pole region, which is caused by SA — > qq/ll Via a heavy gauge boson Z H exchange, 
always appears near ~ mz H /2. The mass of Zh is determined by the parameter / in 
the LRTH model. Changing the value of / would lead to the shift of the pole region, as 
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FIG. 11: Contour plot of the relic density in the wia vs. 82 plane. LRTH relation of 82 = 28\ is 
imposed. The region enclosed by two contour lines corresponds to the WMAP 3 a region. 




FIG. 12: The left plot (a) shows the relic density Qh 2 vs. ma for / = 500 GeV (solid curve), 
600 GeV (dashed curve) , and 700 GeV (dash-dotted curve) . The right plot (b) shows the contour 
plot of the relic density in the m$ vs. / plane. The band in plot (a) and the region enclosed by 
two contour lines in plot (b) are the WMAP 3 a regions. The mass splittings 82 = 28\ = 2 GeV 
is chosen following the LRTH relation. 



evident in Fig. [T2] (a). For larger values of /, the pole region shifts to larger values of rrig. 
Fig. [T2] (b) shows the WMAP 3 a region in the vs. / plane. The region on the left 
which is independent of / corresponds to the high mass bulk region I, while the region on 
the right where rrig grows with / corresponds to the high mass pole region. 

Fig. [13] shows the contour plot of the dark matter relic density in the LRTH model 
parameter space / vs. fi. For fi=0 GeV, where the mass of the dark matter is mainly 
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signal [GeV] 

FIG. 13: Contour plot for the dark matter relic density in the fi vs. / plane for the LRTH 
model. The mass splitting 62 = 28\ = 2 GeV (solid curves) and 10 GeV (dashed curves) are 
chosen following the LRTH relation. The region enclosed by two contour lines corresponds to the 
WMAP 3 a region. 



provided by the CW potential in the LRTH model, we obtain a region in /, / > 1000 GeV, 
that falls into the WMAP 3 a window. This is encouraging that with the minimal setup 
of the LRTH model without introducing an extra mass parameter fi 2 for the dark matter, 
it could accommodate the amount of cold dark matter in the Universe. 



V. CONCLUSION 

The twin Higgs mechanism provides an alternative solution to the little hierarchy 
problem. When implemented in the left-right models, there is a natural candidate for 
WIMP dark matter. One SU(2)l Higgs doublet couples only to the gauge sector while 
it does not couple to the SM fermion sector. The stability of the lightest component in 
this SU(2)i Higgs doublet is protected by a matter parity. The dark matter in the LRTH 
model is similar to that in the IHDM. However, it has its unique features due to the extra 
particles in the model. 

In this paper, we analyzed the dark matter relic density in the LRTH model when the 
mass splittings between S, A and hf follow the approximate relation 5 2 = 25i. We also 
generalize our results to the case of the IHDM when such a mass splitting relation is not 
imposed. We found that there are two distinctive mass regions for m§ where the relic 
density falls into the WMAP window: (A) low mass region, and (B) high mass region. For 
the low mass region, there are coannihilation pole regions when at least one of 5\ or 5 2 is 
small. There is also a bulk WW region around m§ ~ 70 GeV for small values of 62 and 
a bulk bb region for larger values of 62 ■ Since the relevant particle masses and interactions 
that entered the relic density calculation also appear in the IHDM, our results also apply 



16 



to the dark matter relic density in the IHDM. 

In the high mass region, there are high mass bulk regions (I) and (II) due to SS — > 
WW/ZZ for 5\ t 2 ~ 15 GeV. There is also a high mass pole region due to the coannihilation 
of S and A via Zh exchange. While the pole region only appears in LRTH, the bulk region 
also exists in the general IHDM. 

In our analysis, we have assumed that the scalar component S of h® is lighter than the 
pseudoscalar component A. Similar numerical results can be obtained for a pseudoscalar 
dark matter candidate. 

In summary, there are relatively large regions of parameter space in the LRTH model 
(as well as in the IHDM) that could accommodate the amount of cold dark matter in 
the Universe. Even in the minimal setup of the LRTH model with no additional mass 
parameter introduced for the dark matter, the relic density is still sufficient to close the 
Universe. The direct and indirect detection for such a dark matter candidate at current 



and future dark matter detection experiments can be found in Ref. [15 
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